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A tapered optical nanofiber simultaneously used to trap and optically interface of cold atoms 
through evanescent fields constitutes a new and well controllable hybrid quantum system. The 
atoms are trapped in two parallel f D optical lattices generated by suitable far blue and red detuned 
evanescent field modes very close to opposite sides of the nanofiber surface. Collective electronic 
excitations (excitons) of each of the optical lattices are resonantly coupled to the second lattice 
forming symmetric and antisymmetric common excitons. In contrast to the inverse cube dependence 
of the individual atomic dipole-dipole interaction, we analytically find an exponentially decaying 
coupling strength with distance between the lattices. The resulting symmetric (bright) excitons 
strongly interact with the resonant nanofiber photons to form fiber polaritons, which can be observed 
through linear optical spectra. For large enough wave vectors the polariton decay rate to free space 
is strongly reduced, which should render this system ideal for the realization of long range quantum 
communication between atomic ensembles. 

PACS numbers: 37.10.Jk, 42.50.-p, 71.35.-y 
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I. INTRODUCTION 

Effects related to light-matter interactions continue to 
be of big interests in the main present researches, due 
to their importance for fundamental physics and appli- 
cations. Manipulation and trapping of atoms are few of 
the sequences of such interactions [l| , and which serve as 
a framework for the development of quantum informa- 
tion technologies, e.g. the implementation of quantum 
memories, optical communications, and quantum compu- 
tations Furthermore, trapping and guiding of atoms 
are useful tools for atom optics and atom interferometry 
Q . Historically the system that take must attention was 
cavity quantum electrodynamics, which offers a power- 
ful tool for controlling atom-photon interactions and 
even achieving the trapping of a single atom in a cavity 

@ . Recently, a Bose-Einstein condensate of an ultra- 
cold atom gas is fabricated between optical cavity mir- 
rors, and were the strong electric dipole coupling regime 
achieved 0,11. 

Different directions of researches in the past few years 
have been opened for the light-matter coupling which rest 
on optical fibers. In hollow core fiber the atoms are fun- 
neled into a capillary in the center of the fiber and they 
couple to the guided fiber photons, where this s;Y^tem can 
guide and confine both atoms and photons |fll - [l3| . Other 
setup which is of importance for our present work is ta- 
pered optical fibers with a nanofiber waist, where here 
the atoms are outside the fiber and couple to the evanes- 
cent field surrounding the fiber [11 This system 
combining both atomic and solid state devices and can 
be considered as hybrid quantum system, where trapping 
and optically interfacing of atoms can be simultaneously 
achieved. Moreover, the strong evanescent field around 
the nanofiber is efficiently used for trapping, detecting 
and manipulating of a single atom [l^ . 

The guided modes of ultrathin optical fibers, with di- 



ameters smaller than the wavelength of the guided light, 
exhibit strong transverse confinement and pronounced 
evanescent field. Interference of two color evanescent 
fields surrounding an optical nanofiber give rise to an 
array of optical microtraps [l7| . The system realized 
recently for cesium atoms interacting with a multicolor 
evanescent field surrounding an optical nanofiber [l8| . 
where the atoms are localized in a one dimensional op- 
tical lattice parallel to the nanofiber, which shown to be 
efficiently interrogated with a resonant light field sent 
through the nanofiber. Moreover, conventional optical 
lattices are formed of counter propagating laser beams 
to get standing wave in free space, and then the ultra- 
cold atoms are loaded [l^ . This system well described by 
Bose-Hubbard model which predicts the quantum phase 
transition from the superfluid into the Mott insulator 
phase with a fixed number of atoms per site 20] . 

In the present paper we investigate a simplified model 
related to the setup presented in [18j. The hybrid quan- 
tum system build of a tapered optical nanofiber and two 
optical lattices. The one dimensional optical lattices 
formed by off resonance nanofiber modes appear paral- 
lel to each other localized at two opposite sides of the 
nanofiber. We consider the case of a single two-level atom 
per lattice site. The system is excited and probed by pho- 
tons in a fiber mode resonant with the atomic transition. 
Due to resonant dipole-dipole interactions the electronic 
excitation dynamics for each chain is dominated by the 
formation of collective electronic excitations (excitons) 
incorporating the optical lattice translational symmetry 
|2l| - |23{ . Naturally excitons in the two chains interact di- 
rectly and via the resonant fiber mode. In the limit of 
long wave length excitons and for a distance of the two 
chains larger than the lattice constant, the interaction 
among the two optical lattices can be calculated explicitly 
and strongly deviates from single particle dipole-dipole 
coupling. Eventually interlattice interactions give rise to 
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dark and bright excitons which mix with the nanofiber 
photons in the strong couphng regime to form fiber po- 
laritons [23 |. 

The paper is organized as follows. In section 2 we 
derive the exciton dispersions in two parallel optical lat- 
tices, and the exciton damping rate is presented. The 
excitons coupling to nanofiber photons is obtained in sec- 
tion 3. The strong coupling regime is treated in section 
4 to give fiber polaritons. The transmission spectrum is 
calculated in section 5, and a summary appears in sec- 
tion 6. The appendix includes rigorous calculations of 
the itnra and inter lattices exciton dynamical matrices. 



II. EXCITONS IN TWO PARALLEL ID 
OPTICAL LATTICES 

We consider two identical parallel one-dimensional op- 
tical lattices with one atom per site, where the lattice 
constant is a and they are separated by a distance d. 
The optical lattice is located at a distance b from the 
fiber surface, where the two optical lattices are at two 
opposite sides of the fiber, as seen in figure (1), which is 
related to the setup in jl^ . The total Hamiltonian of the 
hybrid quantum system is H = H^x + Hph + Hin , where 
Hex is the electronic excitations Hamiltonian in the two 
optical lattices, Hph is the fiber photons Hamiltonian, 
and Hin is the coupling Hamiltonian between the fiber 
photons and the electronic excitations in the two optical 
lattices. First, we concentrate in the dynamics of elec- 
tronic excitations in the two optical lattices. We examine 
the formation of hybrid collective electronic excitations 
(excitons) among the two lattices. An electronic excita- 
tion delocalizes in a lattice due to resonant dipole-dipole 
interactions and form an exciton which is a coherent state 
with wave number k j2I, 22] . Then we check if such an 
exciton can transfer between the two lattices due to res- 
onant dipole-dipole interactions, and if it can coherently 
oscillate between the two lattices to form a collective elec- 
tronic excitation in the whole system of the two optical 
lattices. 
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FIG. 1: Two One-dimensional optical lattices with one atom 
per site, which are located parallel to a fiber at two opposite 
sides. An incident field is sent from the left side and the 
transmitted field is observed at the right side. The lattice 
constant is a and the distance between the lattices is d. Here 
we assume d ^ a. The transition dipole ^ which makes an 
angle 9 with the lattice direction is also seen. 



The electronic excitation Hamiltonian in the two opti- 
cal lattices is 

Hex = '^Ea B\^Bna + ^ Jnm B^^Braf}- (1) 
n.OL nm,al3 

Here (n, m) run over all the lattice sites, and (a, /3) stands 
for the two lattices (1, 2). The atomic transition energy is 
Ea — fvjjA for two level atoms, which is the same at the 
two lattices. iJ^Q,, Bna are the creation and annihilation 
operators of an electronic excitation at site n in lattice a, 
respectively. We consider the case of a single electronic 
excitation or low excitation density, then the operators 
obey bosonic commutation relations. The coupling pa- 
rameter is for resonant dipole-dipole interactions, 
and give rise to excitation transfer among two sites. The 
transfer term include transfer among atoms in the same 
lattice, and J^^j ^'^d transfer among atoms at dif- 
ferent lattices, J}^ and J"^- 

The Hamiltonian can be diagonalized in the lattice 
sites in using the transformation 



B,, 



(2) 



where N is the number of lattice sites, which is 
taken to be large number, and k is the wave number 
which takes the discrete values k = 'WaP^ with [p — 
0, ±1,±2,--- ,±N/2), in using periodic boundary con- 
dition. Here x" is the position of site n in lattice a. 
The new states represent collective electronic excitations 
which are called excitons and they are waves that 
propagate to the left and the right of the lattice with 
wave number k which is a good quantum number. Now 
the Hamiltonian reads 

Hex = E ^L^kc + J2 J'^^i^) ^L^kp, (3) 

k,a k,af3 

which is diagonal in k. We defined the exciton dynamical 
matrix by 



(4) 



and we used = J°'^{L), with L = x^^ ~ x", where 
the dipole-dipole interaction is a function of the distance 
between the two atoms. 

Next we diagonalize the above Hamiltonian relative to 
the two lattice indexes, (a,/?), by applying the transfor- 
mation 



Bku — 



Bki ± 



fc2 



V2 



(5) 



The new states represent entangled states between ex- 
citons from the two lattices, which are symmetric and 
antisymmetric states. The symmetric state is denoted 
by {v = s) and takes the plus sign the antisymmet- 
ric state is denoted by {v = a) and takes the minus sign 
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(— ). The Hamiltonian casts into the diagonal form 

H,, = Y,KAk)BlBk,., (6) 



where the eigenenergies are 



El^{k) = EA + J{k) + J'{k), 
El^{k) = Ea + J{k) - J'{k). 



(7) 



For identical lattices, inside the same lattice we have 

J{k) = J^\k) = J^^{k), (8) 

and among different lattices we have 

j'{k) = J^^ik) = J^Hk). (9) 

Our main task now is to calculate the two exciton dy- 
namical matrices, J{k) and J'{k). 



A. Exciton Dispersion 

The resonant dipole-dipole interaction between two 
atoms which arc separated by a distance R, and of tran- 
sition dipole is defined by 



J{R) = 



1 



3(/i • R)2 



47reo 



|R|5 



(10) 



Such interaction holds in the limit of |R| < \a, where 
Ea = hc/\A, otherwise one needs to include radiative 
corrections. For the dipole moment wc take the general 
case of /2 = {^xi fJ-yT fJ'z)- The two optical lattices are par- 
allel to the X axis, one with zi = i) and sites at (la, 0,0), 
and the other lattice with Z2 = d and sites at {la,0,d), 
where I is an integer that runs over the whole lattice sites. 
See figure (1). 

In appendix A.l we calculate the exciton dynamical 
matrix for interactions inside the same lattice, J{k). In 
the limit of long wavelength excitons, that is fca <C 1, we 
get 



J{k) 



C(3) My + - 
27reo a'^ 



(11) 



Here the y and z components, which are normal to the 
lattice direction, are repulsive; and the x component, 
which is parallel to the lattice direction, is attractive. 
Note that using the nearest neighbor interaction in this 

limit gives J =^^ii4^. 

Appendix A. 2 includes the calculation for the exciton 
dynamical matrix for interaction among the two optical 
lattices, J' (A;). The result for the limit of long wavelength 
excitons, that is A;a <C 1, and large distance between the 
two lattice kd^ 1, namely d » a, is given by 



^ -fcdj 2_ 2 

arf2 47reo ^ ^ 3 kd 



27T {kd) 



(12) 



Here the x and y terms arc repulsive, as they parallel 
among the different lattices, but as the x component 
is parallel to the lattice direction and the y component 
is normal to it, the contribution of the x component is 
much larger than the y one. The z term is attractive, 
even though these components are normal to the lattice 
direction they are parallel for different lattices. Approx- 
imately the last term is much smaller than the first two, 
(or we can assume fiy = 0), we get 



J'{ka -> 0, fed -)> oo) 



27r(M)3/2 _^ 



ad"^ 47reo 



(13) 

The main resiilt of the present derivation is that the in- 
teraction decay exponentially with the distance between 
the two lattices. In place of the inverse cube dependence 
of the interaction between two dipolcs, the collective ef- 
fect in each lattice gives coupling between two excitons 
which is decay exponentially with the distance. 
The exciton dispersions, {v = a,s), for Hy = 0, are 



E^^ik) ^ Ea + 



± 



C(3) 
2TTeoa^ 
2i {kdfl^ 
ad? 47reo 



{m'-2/x'} 

e-'^'^K-M^}. (14) 



For this case of a dipole in the (x — z) plane, with 
A* = (Mx)0, /U^) = /i(cos^, O,sin0), where 9 is the angle 
between the dipole and the lattice direction, as seen in 
figure (1), we get 

E'^,,{k, 0) ^ Ea + {1-3 cos^ 0} 



± ^M!^!e-^V{2cos^^-l}.(15) 

The intralattice exciton dynamical matrix J{k) is much 
larger than the interlatticcs one J'(k), in the present 
limit. The splitting between the two exciton dispersions 
is 



A{k,0) 



V2^ jkdf''^ 
ad? 47reo 



-kd,,2 



IJ? [2 cos^e- I) , (16) 



which is the coupling parameter for the coherent oscil- 
lation of the exciton between the two lattices. Such os- 
cillation is possible only if the splitting is larger than 
the exciton damping rate. The symmetric-antisymmetric 
splitting vanishes at ^ = 45°, where we get the energy 
^(61 = 45°) K Ea - At 6' w 54.7° the on-lattice 

term vanishes, where 1 — 3 cos^ ^ = 0, and we get 



i?;;(fc,0«54.7°) ^ Eat 



2^ {kdf'n'' 
3ad'^ 47reo 



-kd 



.(17) 



B. Exciton damping rate 

In spite the existence of the nanofiber, the optical lat- 
tices are located in free space, (in the next section we 
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consider the coupling to the fiber photons). Therefore, 
we present here the damping rate of the exciton into 
free space. For a single one dimensional optical lattice 
with one atom per site we calculated in '25'] the radiative 
damping rate into free space. The exciton damping rate 
is given by 



{hckf 



{1 



cos 



(2 cos^ 



sm 



(18) 



which is a function of k and 9. Here the exciton en- 
ergy is Eex{k) = Ea + J{k), which is as the above result 
but in neglecting J'{k). This result needs to be com- 
pared with the single atom the damping rate, which is 

1,2 p3 



As we discussed in more details in (2f 



the exciton damping rate is found to be much larger than 
for small wave numbers, that is fca ^ 1, and the exci- 
tons can be considered as superradiant states. For larger 
wave numbers, the damping rate of excitons with polar- 
izations orthogonal to the lattice increases, while for po- 
larizations parallel to the lattice decreases and becomes 
zero at a critical wave number. As the spontaneously 
emitted photon has wave vector component parallel to 
the lattice equals to the exciton wave number, the crit- 
ical wave number is obtained from the conservation en- 
ergy condition i?|^(fcc) = hck^- Beyond this critical wave 
number, excitons become metastable with zero radiative 
damping rate, and they can propagate in the lattice with- 
out radiative decay. 

In the present system we have two optical lattices, and 
we showed above that their mutual coupling give rise to 
symmetric and antisymmetric excitons. Direct calcula- 
tions yield that the antisymmetric excitons are dark with 
zero damping rate, that is F^^ = 0, while the symmetric 
excitons are bright with twice the single exciton damp- 
ing rate, that is F^^. = 2Fe2;. Note that according to the 
previous discussion, the symmetric exciton damping rate 
also become zero beyond the critical wave number kc- 
The dark excitons can not be excited directly through 
the coupling to the fiber photons and needs other meth- 
ods to excite them. As we explain in the next section, the 
symmetric excitons can be excited directly through the 
electric dipole interaction with the fiber photons, where 
this coupling appears also for excitons with wave number 
beyond the critical wave number. This regime is useful 
for long range transmission through the system, and the 
only source of damping will be the fiber photon decay. 



III. FIBER PHOTON-EXCITON 
INTERACTIONS 

Now we treat the fiber photons and their coupling to 
the two optical lattices. Ri goro us treatment of tapered 
nanofiber modes is given in |14| . but here we use a sim- 
plified picture. The one-dimensional optical fiber modes 



are given by the Hamiltonian 



Hph = y^^hwphjk) 



(19) 



where and are the fiber photon creation and anni- 
hilation operators with wave number k, respectively. The 
photon dispersion is taken to be of the simplified form 



ujphik) = — ^/cj 



(20) 



where the wave number kg results of the fiber transverse 
confinement, and e is the fiber average dielectric constant. 
The electric field operator outside the fiber is 



E(r) = I 



(21) 

where e is the photon linear polarization unit vector, V 
is the normalization volume, and u(r) is the mode func- 
tion which includes the fiber field complexity. The fiber 
photons can decay into free space with a given damping 
rate fixed by the nanofiber quality. 

The atomic transition dipole operator is fl — 
t^J^n.a {Ena + Bj^a)- The matter-light couphng is given 
by the electric dipole interaction Hin — — /i • E, and in 
the rotating wave approximation for linear polarizations, 
we get 



Hi, 



k,n,a 



{n-e)uib) 



^-"}, (22) 



where in using the inverse transformation of equation ([2|) , 
we have 



Hi, 



/ hwph{k)^'^ 



kl.: 



2eQVN 



Bic, e-'('=-')^"}. (23) 



In using J2n ^'^'^ — ^ki, we obtain 



Hin = -2^u(6)^ 



lhwphik)^i^N 



k.a 



2eoV 



(24) 

We represent the exciton operators in terms of symmetric 
and antisymmetric operators by using the inverse of the 
transformation ([S]), to get 



Hin = -i^u(6)^ 



lhuJphik)^l^N 



eoV 



{akBl-alBks} 



(25) 
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It is clear that only the symmetric excitons are coupled 
to the fiber photons, where the antisymmetric ones are 
dark and decouple to the photons. 
In defining the coupling parameter 



^fk = -iu{b)\ 



eoSa 



(26) 



where S is the mode cross area, then we can write 

H,u = 51 {^/^ "''^L + 4Bks } . (27) 
fc 

The photon dispersion and the exciton-photon coupling 
is taken to be 9 independent. 

IV. FIBER POLARITONS 

The total Hamiltonian is given by 

k 

+ ujph{k) a\ak + fk o-kBl^ + Jl aX^ks^ , (28) 

where due to translational symmetry along the lattice 
and fiber axis, the Hamiltonian is separated for each fc. 

In the strong coupling regime, where the coupling is 
larger than both the exciton and photon line width, we 
define polaritons by diagonalizing the Hamiltonian [2l|, 
[23 | , to get the total Hamiltonian by 



(29) 



kr k 

with the polariton dispersions 



where 

A{k,9) = ^S^{k,9) + \fk\^, 
with the detuning 

gf^f^ ^ ^phjk) - ujl^{k,9) 



(30) 



(31) 



(32) 



The polariton operators are a coherent superposition of 
symmetric excitons and photons, where 



At 



X^{k,9) Bks+Y^{k,9) ak, 
and where the amplitudes are given by 



(33) 



X^{k,9) = 
Y^{k,9) = 



'A(fc,6') =F'5(fc,e') 



2A(fc,6' 
fk 



^2A{k,9){A{k,9)T6{k,9)) 



(34) 



Here we present the results for a system with the fol- 
lowing parameters: the lattice constant is a = 1000 A, 
the transition energy is Ea = 1 eV, the transition dipole 
is /i = 1 eA of angle 9 = 90° with the lattice axis. The 
distance between the two optical lattices is taken to be 
d = 10a. The fiber dielectric constant is e = 3, with the 
mode cross area of 5* = 47ra^, and the mode function at 
the lattice position is taken to be u{b) = 0.1. The cavity 
mode energy at fc = is taken to be in resonance with 
the free atom transition energy, that is Eph{k = 0) = i?^, 
then fco = EA^/e/hc. The polariton dispersions are plot- 
ted in figure (2) relative to the free atom transition, and 
the symmetric exciton and photon weights are plotted 
in figure (3). The exciton-photon intersection point is 
taken here to be at fc = 0, where the polaritons split 
by the Rabi splitting, and they are half exciton and half 
photon. For large fc the upper branch becomes photonic, 
and the lower becomes excitonic. 
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FIG. 2: The full lines are the upper and lower polariton angu- 
lar frequencies relative to the free atom transition, uipoi — oja 
vs. wave number fc. The dashed line is the symmetric exciton 
dispersion and the dashed parabola is for the fiber photon 
dispersion. 



V. FIBER TRANSMISSION SPECTRUM 

To get the linear optical spectra f^TIl, we consider an 
incident field from the far left side of the fiber, and we 
calculate the transmission T{k, 9) and reflection i?(fc, 9) 
spectra, with the absorption spectrum A(k, 9). The cou- 
pling of the fiber field to the external field at the two 
far edges of the fiber is included in the parameter 7, and 
the fiber photon damping rate is taken to be Tph which 
is assumed to be fc independent. The symmetric exciton 
damping rate is included here phenomenologically, and 
is taken to be Fg^(fc, 6*) as discussed in section 2. A. The 
polariton damping rate is taken to be 

^toi^k,9) = TlAk,9) \X^{k,9)\^+Tp^ \Y^{kM^- (35) 
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FIG. 3: The excitonic and photonic weights vs. k. For the 
lower polariton branch the dashed Une is for the excitonic 
weight, and the full line for the photonic weight. For the 
upper polariton branch the dashed line is for the photonic 
weight, and the full line for the excitonic weight. 



The linear optical spectra are calculated and given by 

1 



and 



with 



where 



R{uj,k,0) = — 



T{uj,k,( 



+ ij A(a;,fc,^)|2' 

>- |A(:..A-,/?)l2 



|l + i7 A(a;,A;,^)|2' 



R{uj, k, e) + T{io, k, 6) + A{uj, k, 6) = 1, 



A{Lo,k,9)=Y,^ 



Y-ik,e)\'^ 



CO — CO, 



pol 



{k,ey 



(36) 



(37) 



(38) 



(39) 



where u,;jk, 0) = uj;jk, 0) - iVfjk, 6). 

In the hncar spectra plots we use for the photon a 
small damping rate of hVph = 10~^° eV, and the fiber 
edge coupling parameter is taken to have a bandwidth 
of — 10^^ eV. In figure (4) we plot the transmission 
spectra for a wave number of A; = 10~^ where we 

get two transmission peaks correspond to the two polari- 
ton branches which are separated by the Rabi splitting. 
Now we take an incident field of a wide bandwidth of 
h'y = 10~^ eV. In figure (5) we plot the transmission 
spectra for a wave number of A; = 10~^ where here 
the band gab between the two peaks becomes a dip in 
the transmission spectrum. Figure (6) is a close look 
at the minimum which show a small shift relative to the 
free atom transition, and that results of the dipole-dipole 
interactions. Here the symmetric exciton line width is 
hri^ = 2.32 X 10"* eV, which is larger than the free 
atom line width of UTa = 2.5 x 10~^ eV. In figure 
(7) we plot the transmission spectra for wave number of 
k = 10^5 A-\ and in figures (8-9) for A: = 5 x 10"^ 



where in figure (8) the plot is around the upper polariton 
branch, and in figure (9) is around the lower polariton 
branch which is also seen in figure (8) as a small peak 
around zero. It is clear that for large wave numbers at 
the lower branch we get small transmission, where for 
large k the excitons and photons with the same wave 
number have large detuning, namely they are off reso- 
nance. At the upper branch we get large transmission. 
As we mentioned previously, the upper branch becomes 
photonic and the lower one excitonic. 
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FIG. 4: The transmission spectrum, for fe = 10 ^ A ^, and 
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FIG. 5: The transmission spectrum, for fc = 10 ^ A ^, and 
= 10-* eV. 



VI. SUMMARY 

We studied a recently demonstrated integrated hy- 
brid atom-photon quantum system based on a tapered 
nanofiber and cold atoms. The nanofiber is used to im- 
plement two arrays of nanotraps for the atoms using off- 
resonant photons as well as to spatially confine resonant 
photons efficient to ensure strong single photon coupling 
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FIG. 9: The transmission spectrum, for fc = 5xl0 ''A ^, 
around the lower polariton branch. 




FIG. 7: The transmission spectrum, for = 10 ^ A 



0.6 

E 

CO 

§0.4 

































1 


2ru 1 3 4 



FIG. 8: The transmission spectrum, for = 5 x 10 ^ A 
around the upper polariton branch. 



herently mix with the nanofiber photons to form fiber 
polaritons, which can be seen in the transmission spec- 
trum of the nanofiber. 

The present system allow one to investigate and to 
deeply understand one dimensional collective electronic 
excitations (excitons), where it is possible to excite, store, 
and readout information in the form of coherent elec- 
tronic excitations. As was shown in our previous work 
[2^, the damping rate of small wave number excitons 
in one dimensional lattice is much larger than a single 
atom excitation, and decay fast into free space. But be- 
yond a critical wave number it is shown that one dimen- 
sional excitons became metastable with infinite radiative 
life time. In this regime, excitons can not be excited 
by free space photons, but can be easily excited by cou- 
pling to fiber photons. These excitons can serve for infor- 
mation storage devices and for communications through 
energy transfer. In the present paper the discussion is 
limited for the linear regime. Nonlinear processes in- 
duced by electronic excitation saturation effects give rise 
for more interesting physics, which is our next task. At 
this point we have neglected the interaction of the atoms 
via light forces induced by multiple scattering of pho- 
tons [26[, which can be significantly enhanced through 
the transverse mode confinement and the fact that the 
atoms form Bragg like structures. This could lead to col- 
lective oscillations and instabilities of the lattices but also 
can be used to engineer long distance motional couplings 
between ensembles. 



to the two atomic ensembles. Interestingly the effective 
coupling parameter of excitons in the two chains decays 
exponentially with their distance but it is still sufficient 
to split the energies of symmetric and antisymmetric ex- 
citons. Only the symmetric excitons are bright and co- 



Acknowledgments 



The work was supported by the Austrian Science 
Funds (FWF) via the project (P21101), and the Eranet 
NanoSci-E-t- project NOI (I269-N16). 



8 



Appendix A: Exciton dynamical matrix calculations 

In this appendix we calculate the exciton dynamical 
matrix for the two cases of intralattice interactions, J(fc), 
and interactions among interlattices, J'(fc). 



1. Intralattice interactions 

For interactions inside the same lattice, from equation 
T0| . we have the resonant dipole-dipole interaction 



J{1) 



1 + 



47reo a'^P 
and in momentum space reads 



(Al) 



J(fc) 



1 ^4+^4-2(4"-^' e^'^^^ 



47reo o,^ 

1— — 00 

1 + Mz ^ 2/1^ 2 cos{kal) 



47ren ^ P 

' 1=1 



(A2) 



the prime indicate that I runs from — oo to +oo in ex- 
cluding I = which is for self interactions. In the limit 
of long wavelength excitons, or small wave number, that 
is fca <C 1, we get 



1=1 



(A3) 



27reo ■ ■ a 
where the Riemann zeta function gives C(3) ~ 1.202. 

2. Interlattices interactions 

In order to calculate the exciton dynamical matrix for 
interactions among the two lattices, we use the Ewald's 
method [131, which is applied here for excitons in one 
dimensional optical lattices. 

The resonant dipole-dipole interaction reads 



(A4) 



where = x, y, z). The calculations give 
1 3a^P 



(a2/2+rf2)3/2 (^2^2 _^ ^2)5/2 ' 
1 



(a2/2 +^2)3/2' 
1 



' {aH^ + d2)3/2 (^2^2 + ^2)5/2 ' 

) = D (l)^ — 

' ""^^ (a2p+d2)5/2' 

) = Dy,{l) = Dy.il) = D.yil) = 0. (A5) 



We need to calculate the interactions in momentum space 

(A6) 



where 



+00 



l= — oc 



We define the function 



S{k) ^ E 



l=-c 



(02^2+^2)5/2' 



(A7) 



(A8) 



then the matrix elements of the exciton dynamical matrix 
can be derived from it by 



92 



Dyy{k) - ( + ) ^W, 



a/c2 

D,,{k) = DlJk) = t3d^S{k). (A9) 

ok 

The function S{k) can be written in a form of an inte- 
gral with exponential decay, in place of the l/R^ depen- 
dence, in using the identity 



3\/7r Jo 
Then we arrive at 



4 , I 



,5/2 " 



(AlO) 



S{k)^-^ / dtt^'^e-""' y 6-'^'''*+^'='". (All) 



In using the relation 

+ 00 , +00 



^ a^/t ^ ^ 



n— — oo 



We convert one summation into another. The first in- 
cludes oscillations, and then decays slowly for large sum- 
mation index I. While the second decays faster with in- 
creasing the summation index n. Such a change in the 
series found to be useful for long wavelength excitons, as 
we show in the following. In doing this, we obtain 

S{k) = —Y / i e-'*'*e-^H+'^)'. (A13) 

The integral gives 

/ dt t ^ 2^K2{2y//3^), (A14) 

Jo 7 
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where Kn (x) is the modified Bessel function of the second 
kind of order n. Finally we have 



4 

3a 



+ 00 n ( 



'2d ( 



(A15) 

We interest in the limit of long wavelength excitons, 
that is /ca ^ 1. In this limit the dominant contribution 
comes from the term with (n = 0), and the other terms 
decay faster, hence we get 



We need to use the relation 
d 1 

— Kn{x) = --[Kn-l{x) + Kn+l{x) 



(A16) 



(A17) 



Here we use the second main approximation, which is for 
large d where fcrf 1, and as fca ^ 1, we are in the limit 
of a. For the function Kn{kd) we use the asymptotic 
expansion 



Kn{kd oo) 



2kd 



-kd 



(A18) 



which has the same behavior for different n. In this limit 
we get the matrix elements 



Dxx{ka 0,kd^ oo) 
Dyy{ka — >■ 0, fed — )• oo) 
Dzz{ka -> 0, fed -> oo) 
Dxz{ka -^Q,kd^ oo) 

Finally we obtain 

J'{ka -^0,kd^ oo) 



ad2 
4 

V2n 



{kdfl'' e-'^^ 



-2i 



{^kdfl^ e 



■kd 

(A19) 



2^ {kdf^ 
ad'^ 47reo 



-kd 



[1] H. J. Mctcalf, and P. van dcr Stratcn, Laser Cooling and 

Trapping, (Springer, NY, 1999). 
[2] D. Bouwmeester, A. K. Ekert, and A. Zeilinger, The 

physics of Quantum Information, (Springer, NY, 2000). 
[3] P. Mcystrc, Atom Optics, (Springer, NY, 2001). 
[4] S. Harochc, J. M. Raimond, Exploring the Quantum: 

Atoms, Cavities and Photons, (Oxford UP, UK, 2006). 
[5] J. Ye, D. W. Vernooy, and H. J. Kimble, Phys. Rev. Lett. 

83, 4987 (1999). 
[6] P. W. H. Pinksc, T. Fischer, P. Maunz, and G. Rempe, 

Nature 404, 365 (2000). 
[7] F. Brcnncckc, T. Donncr, S. Ritter, T. Bourdel, M. Khl, 

and T. Esslinger, Nature 450, 268 (2007). 
[8] Y. Colombe, T. Steinmetz, G. Dubois, F. Linke, D. 

Hunger, and J. Reichel, Nature 450, 272 (2007). 
[9] J. Verdu, H. Zoubi, C. Koller, J. Majer, H. Ritsch, and 

J. Schmiedmayer, Phys. Rev. Lett. 103, 043603 (2009). 
[10] A. Imamoglu, Phys. Rev. Lett. 102, 083602 (2009). 
[11] C. A. Christonscn, S. Will, M. Saba, G. B. Jo, Y. I. Shin, 

W. Kctterlc, and D. Pritchard, Phys. Rev. A 78, 033429 

(2008). 

[12] M. Bajcsy, S. Hofferberth, V. Balic, T. Peyronel, M. 

Hafezi, A. S. Zibrov, V. Vuletic, and M. D. Lukin, Phys. 

Rev. Lett. 102, 203902 (2009). 
[13] A. R. Bhagwat, and A. L. Gaeta, Opt. Express 16, 5035 

(2008). 

[14] G. Sague, E. Vetsch, W. Alt, D. Meschede, and A. 
Rauschenbeutel, Phys. Rev. Lett. 99, 163602 (2007). 



[15] K. P. Nayak, P. N. Melentiev, M. Morinaga, F. L. Kien, 
V. I. Balykin, and K. Hakuta, Opt. Express 15, 5431 

(2007) . 

[16] K. P. Nayak, and K. Hakuta, New J. Phys. 10, 053003 

(2008) . 

[17] G. Saguc, A. Baadc, and A. Rausciicnbcutcl, New J. 

Phys. 10, 113008 (2008). 
[18] E. Vetsch, D. Reitz, G. Sague, R. Schmidt, S. T. 

Dawkins, and A. Rauschenbeutel, Phys. Rev. Lett. 104, 

203603 (2010). 

[19] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[20] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. ZoUer, Phys. Rev. Lett. 81, 3108 (1998). 
[21] H. Zoubi, and H. Ritsch, Phys. Rev. A 76, 013817 (2007). 
[22] H. Zoubi, and H. Ritsch, Europhys. Lett. 82, 14001 

(2008) . 

[23] H. Zoubi, and H. Ritsch, New J. Phys. 10, 23001 (2008). 
[24] H. Zoubi, and H. Ritsch, Europhys. Lett. 87, 23001 

(2009) . 

[25] H. Zoubi, and H. Ritsch, Europhys. Lett. 90, 23001 

(2010) . 

[26] J. K. Asboth, H. Ritsch, and P. Domokos, Phys. Rev. 

Lett. 98, 203008 (2007). 
[27] M. Born, and K. Huang, Dynamical Theory of Crystal 

Lattices, (Clarendon Press, Oxford, 1954). 



